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ABSTRACT 

We present the far-infrared (IR) maps of a bipolar planetary nebula (PN), NGC 650, at 24, 70, and 160/im 
taken with the Multiband Imaging Photometer for Spitzer (MIPS) on-board the Spitzer Space Telescope. While 
the two-peak emission structure seen in all MIPS bands suggests the presence of a near edge-on dusty torus, 
the distinct emission structure between the 24/im map and the 70/160/im maps indicates the presence of two 
distinct emission components in the central torus. Based on the spatial correlation of these two far-IR emission 
components with respect to various optical line emission, we conclude that the 24//m emission is largely due to 
the [O IV] line at 25.9/im arising from highly ionized regions behind the ionization front, whereas the 70 and 
160/im emission is due to dust continuum arising from low-temperature dust in the remnant asymptotic giant 
branch (AGB) wind shell. The far-IR nebula structure also suggests that the enhancement of mass loss at the 
end of the AGB phase has occurred isotropically, but has ensued only in the equatorial directions while ceasing 
in the polar directions. The present data also show evidence for the prolate spheroidal distribution of matter 
in this bipolar PN. The AGB mass loss history reconstructed in this PN is thus consistent with what has been 
previously proposed based on the past optical and mid-IR imaging surveys of the post- AGB shells. 
Subject headings: circumstellar matter — infrared: stars — planetary nebulae: individual (NGC 650) — stars: 
mass loss 



1. INTRODUCTION 

NGC 650 (PK 130-10.°I,M 76, L ittle Dumbbell Nebula) is 
a large (~ 300^^: iBalick et alJ ll992) bipolar planetary nebula 
(PN) of the "late butterfly" type (BaHck 1987). The nebula 
structure in the optical consists of the bright rectangular core 
of 95" X 40" (the long side perpendicular to the bipolar axis) 
and a pair of fainter lobes extending ^ 90" and ^ 150" from 
the central star that is attached to the long side of the rectan- 
gular core. The central core was long suspected to be a nearly 
edge-on torus, and the most recent kinematical study has elo- 
quently demonstrated that the core is an inclined torus and the 
lobes ar e blown-bubbles expanding into the polar directions 
jBrvce et al..l996.) . 

Under the framework of the widely-acc epted general inter- 
acting stellar wind (GISW) model (e.g., lBahcklfT987 t). the 
formation of bipolar PNs is understood as a two-step pro- 
cess. First, the progenitor star has to lose its envelope ma- 
terial via mass loss (M ^ 10"^ Mq, v ^ 10 km s"') during the 
asymptotic giant branch (AGB) phase of the evolution. This 
AGB wind leads up to the so-called superwind (M ^ 10""* 
Mq, V 20 km s"') at the end of the AGB phase prior to 
the exhaustion of the envelope material. During the super- 
wind phase, mass loss is expected to cause concentration of 
the ejected matter into the equatorial plane. The AGB mass 
loss therefore results in the equatorially-enhanced circumstel- 
lar shell. 

Then, a hot and tenuous fast wind (M ~ 10"^ Mq, v ^ 10^ 
km s"' ) begins to emanate from the progenitor just prior to 
the beginning of the PN phase when the circumstellar mate- 
rial starts to ionize. This fast wind pushes into the surround- 
ing equatorially-concentrated envelope, and the flow of the 
wind is thus channeled into the polar directions. The prefer- 
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ential wind flow towards the polar directions creates the typi- 
cal bipolar lobes as wind-blown bubbles. The GISW scheme 
is thus capable of producing various PN shapes (both bipolar 
and spheroidal shells) depending primarily on the degree of 
the equatorial enhancement in the surrounding AGB shell. 

Hence, in order for the GISW scheme to work, the 
equatorially-enhanced AGB shell must be present when a fast 
wind is initiated. The GISW model itself, however, does not 
address how to generate the equatorial enhancement in the 
AGB wind shell. The presence of such circumstellar shells 
with a built-in equatorial density enhancement has been con- 
firmed by imaging surveys of proto-PNs in thermal IR dust 
emission from the innermost torus (e.g.. Skinner et al. 1994 
iMeixner et al. 1997; Dayal et al. 1998; Meixneretal. 199^ 
lUeta e t al.ii2001) and dust-scattered star light (e.g., Ueta et a^ 
2000). Based on the results fr om the largest mid -IR and op- 
tical imaging surveys to date (Meixner et al. 1999; Uet a et all 
2000), a so-called "layered shell model" has been proposed to 
explain bot h the bipolar an d elliptical proto-PN morphologies 
(iUeta.2002tlUeta & M'eixner^2003I In this model, the proto- 
PN shells are thought to consist of three generic layers, each 
of which possesses a specific structure reflecting the geome- 
try of mass loss at the time. The outer spherically symmet- 
ric layer represents isotropic mass loss during the early AGB 
phase, while the inner toroidal layer embodies equatorially- 
enhanced mass loss during the superwind epoch at the end of 
the AGB phase. Between these layers there is an intermedi- 
ate spheroidal layer, which results from a gradual transforma- 
tion of the mass loss geometry from isotropic to equatorially- 
enhanced over the course of the AGB phase. 

The unique strength of this layered shell model is its ver- 
satility. Radiative transfer calculations using this model have 
successfully reproduced both bipolar and elliptical proto-PN 
morphologies by simply varying the degree of the equato- 
rial density en hancement (i.e., optica l depth) of the shell 
jMeixner et al J 12002; Uet a & Meixned 12003). The impHca- 
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tion of the layered shell model is that there is no intrin- 
sic distinction between the elliptical and bipolar proto-PN 
shells other than the degree of the equatorial density en- 
hancement. Bipolar PNs would emerge from equatorially- 
enhanced proto-PNs and elliptical PNs would descend from 
proto-PNs of rather isotropic density distribution. The pres- 
ence of the spheroidal intermediate layer is the key to account 
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It is clear that the equatorially-enhanced superwind initiates 
the subsequent aspherical structure development that eventu- 
ally results in the observed complex PN structures. How- 
ever, it is still not yet evident how the shell structure devel- 
opment ensues in the AGB and post-AGB winds. How does 
the spheroidal density distribution arise prior to the equatorial 
enhancement at the end of the AGB phase? Is the post-AGB 
mass loss significant in shell shaping? Does a fast wind sim- 
ply "snow-plow" the surrounding material, tracing the pre- 
existing toroidal density distribution, or contribute to an ad- 
ditional equatorial enhancement? These are only a few ques- 
tions concerning the shell structure development in the cir- 
cumstellar shells of evolved stars. 

Further to understand the role of mass loss in the shell struc- 
ture development during the AGB phase and beyond, the rem- 
nant AGB shells need to be investigated in PNs via sensitive 
far-IR observations of thermal dust emission, since such rem- 
nant AGB shells are very much dispersed and cold. In the fol- 
lowing, we present the results of Spitzer Space Telescope far- 
IR mapping observations of a bipolar PN, NGC 650, to probe 
the history of AGB and post-AGB mass loss in this object. 
Below, we describe the data set and reduction procedure (§2), 
discuss the results (§3), and summarize conclusions (§4). 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 650 was observed with the Multiband Imaging 
Photometer for Spitzer (MIPS; Ri eke et alJ 12004 on-board 
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September 18 as part of GTO first year observations (Pro- 
gram ID 77; AOR KEY 9548032). The MIPS observations 
were performed in the scan-map mode using the medium scan 
rate. The scan was done along four scan legs of 0.°5 long each 
with the cross-scan stepping of 148". The whole scan se- 
quence was repeated twice. The entire three-band dataset was 
retrieved from the Spitzer Science Archive using the Spitzer- 
Pride software tools (Spot/Leopard). 

Data reduction was done mostly with the 2005 September 
30 version of the Mosaicker software provided by the Spitzer 
Science Center (SSC)^ and the Ge Reprocessing Tools (GeRT; 
Ver 041506 of S14 processing)-' . The 24/im map was con- 
structed using the Mosaicker from the Basic Calibrated Data 
(BCD) products that were produced from the raw data via the 
SSC pipeline. Upon making a mosaic map, anomalous pixels 
were removed by using the dual outlier detection algorithm in 
the Mosaicker. 

Since the post-BCD mosaicked maps at 70 and 160/im gen- 
erated by the SSC pipeline showed obvious artifacts, we gen- 
erated these maps directly from the raw data using the GeRT 
and the Mosaicker. First, we corrected bad stimulator re- 
sponse calibration affected by the bright objects (i.e., the tar- 
get), in order to optimally remove the time-dependent respon- 
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sivity variation of the detector This is done by ignoring the 
stimulator flash frames on and near the target. Second, we 
removed the "streakings" due to residual slow response varia- 
tions of the detector by applying a high-pass time filter Since 
we are interested in the extended far-IR structure of the object, 
we did not perform column filtering for the 70/im band data. 
Instead, we used a relatively small median-count value of six 
when applying a high-pass time filter Upon applying the 
time filter, we masked out the target to avoid extended faint 
emission from being mistakenly filtered out. These GeRT 
processes resulted in cleaner BCDs. Then, we constructed 
mosaic maps from the GeRT-processed BCDs using the Mo- 
saicker 

The resulting maps are in 1 ."28 pix"' , 4."92 pix"' , and 7."99 
pix"' for the 24, 70, and 160/im bands, respectively. After 
making artifact-removed maps for all three bands, we sub- 
tracted sky emission, which was estimated at off-source posi- 
tions that are free from background sources. The sky emission 
was determined to be 20, 9, and 27 MJy sr"' respectively in 
the 24, 70, and 160/im bands. These values are consistent 
with the estimated values of the date of observations obtained 
using the Spot software. In order to facilitate direct compari- 
son of the nebula structure between two different MIPS bands, 
the 24 and 70/im images were convolved with the beam of the 
longer wavelength bands and their pixel scales were made to 
match with that of the longer wavelength band. Flux cali- 
bration of the maps were done using the conversion factors 
provided by the SSC. 

3. RESULTS AND DISCUSSION 

3.1. T^GC 650 in the Far-IR 

Figure presents the MIPS mosaicked maps of NGC 650 
at lAjiva (top left), 70/im (middle left), and 160/im (bottom 
left) in grayscale and contours. Also displayed for compari- 
son are optical narrow -band images in the [O III] band at 5007 
A(top right) and HaH-[N II] band at 6600 A(middle and bot- 
tom right) in grayscale, obtained from the Instituto de As- 
trofisica de Canarias (IAC) Morphological Catalog of North- 
ern Galactic Planetary Nebulae Archive."* These optical im- 
ages are overlaid with the far-IR contours. The figure caption 
gives the complete detail of the images. 

In all bands, the nebula core is resolved into two emission 
peaks that delineate the bipolar waist of the optically bright 
core (i.e., the short sides of the core). The southwestern peak 
is brighter than the northeastern peak in all bands. The inten- 
sity of the brighter (southwestern) peak is 72, 63, and 29 MJy 
sr"' while that of the dimmer (northeastern) peak is 61, 35, 
and 23 MJy sr both respectively in the 24, 70, and 160//m 
bands. Previous kinematic studies of the nebula in the optical 
emission lines have shown that the bright core is an expand- 
ing torus (dubbed a "napkin ring") inclined at about 75° with 
respect to the line of sight (e.g., Brvce et al. 1996). A similar 
two-peak emission structure at mid-IR has been commonly 
observed in optically thin, equatorially-enha nced proto-PN 
dust shells that are oriented near edge-on (e.g., Meixner et alj 
1993; Skinner et al. 1994; Uetaet al. 2001; Gledhill & YateJ 
2003). The resolved far-IR structure in this PN is consistent 
with the mid-IR structure of progenitor objects, suggesting 
that the innermost torus of warm dust in a proto-PN shell 
would keep its toroidal shape into the PN phase while the dust 
temperature decreases upon expansion. The two-peaked far- 
IR structure thus represents the limb-brightened edges of a 
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RA offset (arcsec) 

Fig. 1 . — Grayscale MIPS mosaicked maps of NGC 650 at 24fim (top left), lOfim (middle left), and 160/jm (bottom left) overlaid with contours representing 
90 to 10% of the peak (with the 20% interval; in white) and 5, 1, 0.5, and 0.15% of the peak (in black). Also presented are grayscale lAC optical maps in the 
[O III] 5007 Aband (top right) and in the Hq:+[N II] 6600 Aband (middle and bottom right), overlaid with the same far-IR contours as left panels (except for the 
lowest in the 24/im band for clarity). The images are in nominal orientation (N is up, E to the left) and centered at the position of the central star (RA [2000] = 
l''42'"19?656, Dec [2000] = 51°34'32"79; Koomneef & Pottasch 1998). The angular offsets in arcsec are shown by the tickmarks. A circle at the lower right 
comer in the left panels represents the beam size (FWHM) of each band (6", 18", and 40" for the 24, 70, and 160/xm bands, respectively). The number shown 
at the lower left comer in the left panels is the peak intensity in MJy sr"'. The lowest contour is the 3 a level (0.1, 0.9, 3 MJy sr"' respectively for the 24, 70, and 
160/xm bands). 
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near edge-on optically thin dusty torus, at which the column 
density of the far-IR emitting matter is the largest. 

Far-IR emission in the nebula is more extended than the op- 
tically bright core region. Emission is detected at least 10% 
of the peak intensity level in the regions interior to the opti- 
cal inner lobes 90" from the center) in all bands. Fainter 
emission appears to be more extended along the bipolar axis 
beyond the optical inner lobes at shorter wavelengths. In the 
24/im band, emission at ^ 0.5% of the peak level lOa 
detection) is observed in regions corresponding to the outer 
lobes. There is a detached "cap" detected in the optical at the 
tip of t he southeast lobe ( marked in the bottom right panel of 
Fig.ni. lHoraeraD(l200l have reported detection of emission 
from this cap in the Spitzer/Infrared Array Camera (IRAC) 
bands (and even from the opposing one in the northwest lobe). 
Although the 24/im map shows very low-level (^ three a) 
emission in the regions corresponding to the southeast cap, it 
is uncertain whether we have positive detection of the cap or 
not. There appears to be an isolated patch of emission near the 
northwest corner of the panel at (-140", 120"). This is, how- 
ever, not the opposing cap: it is simply a background source. 

Emission of intermediate strengths {r-^ 1 to 20% of the peak) 
shows an interesting morphology of the nebula. In the optical 
image, there is a region of enhanced brightness in both of the 
inner lobes. The bright region in the optical northwest lobe 
seems to be connected to the northeast end of the waist, while 
the bright region in the southeast lobe appears to be attached 
to the southwest end of the waist. On one hand, a similar trend 
is seen in the far-IR maps. The northeast (southwest) peak is 
more elongated towards the northwest (southeast) lobe, corre- 
sponding to the optically bright region. On the other hand, the 
actual location of the northeast (southwest) peak appears to be 
displaced towards the southeast (northwest) direction, which 
is opposite to the direction of the peak elongation mentioned 
above. Figure|2]is presented to show this morphological trend 
by surface brightness profiles. 

While the overall boxy shape of the optically bright core is 
oriented perpendicular to the bipolar axis, the emission peaks 
do not seem to be aligned with the equatorial plane. The elon- 
gation of the peaks is best seen in the 24/im map, while the 
peak displacement is well displayed by the 70/im map. This 
morphological trend is consistent with that ha s been seen in 
the op tical (Balick 1987: Brvce et ^ al.. 1996I) . H2 JKastner et al.1 
119961) . and IRAC Mora et al.l2004l) bands, and gives the neb- 
ula point-symmetric appearance. 

Aperture photometry above the three a level (the lowest 
contour in Figure^ yields flux measurements of 4.5 1 ± 0.04, 
6.04 ±0.30, and 4.83 ±0.95 Jy in the 24, 70, and 160/im 
band, respectively. These values are listed in Table^together 
with the past measurements taken by Infrared Astronomical 
Satellite (IRAS). The measured flux at the 24/im MIPS band 
is about 50% higher than that at the 25/im IRAS band. The 
color correction on the 25/im IRAS flux does not seem to ac- 
count for the entire 50% discrepancy. Since the entire 24/im 
MIPS band is overlapped with the 25/im IRAS band, the dis- 
crepancy cannot be attributed to the presence of strong IR 
fine-structure lines in the 24/im MIPS band that are not cov- 
ered by the 25/im IRAS band. Sensitivity of the 24/im MIPS 
band is superior to that of the 25/tm IRAS band, and the 24/im 
MIPS map shows a very faint extended nebulosity. Thus, the 
discrepancy could be attributed to the IRAS measurements 
missing the faint emission. Aperture photometry in the 24/im 
MIPS band using the estimated measurement threshold at the 
25/im IRAS band yields 3.67 Jy. Considering also the beam 




-6a 6 ba 

RjadiaJ Qflact (arcane ^ 

Fig. 2. — Surface brightness profiles of the emission peaks in the direc- 
tion of the polar axis. The top (bottom) panel is for the northeast (southwest) 
peak. The cuts are made along lines of 25" width parallel to the bipolar axis 
(of PA = -45°). The 24, 70, and 160/im band data are shown respectively by 
dotted, dashed, and solid lines. Positive offsets corresponds to the northwest 
direction as indicated in the panels. While the peak location is shifted to- 
wards southeast (northwest) in the northeast (southwest) peak, the peak itself 
is more elongated towards the opposite, northwest (southeast) direction. 



TABLE 1 

Far-IR Flux Measurements for NGC 650 



Measurements 


Flux 

(Jy) 


Flux 

(Jy) 


Flux 

(Jy) 




24^tm 


70^m 


1 60/im 


MIPS 


4.51 ±0.04 


6.04 ±0.30 


4.83 ±0.95 




25/tm 


60/im 


lOO/im 


IRAS PSC" 
IRAS FSC'' 


2.79 ±0.20 
3.14±0.16 


6.80 ± 1.43 
7.28 ±0.66 


9.29±0.51 
10.3 ±0.9 



"IRAS Point Source Catalog data for IRAS 01391+5119 
''IRAS Faint Source Catalog data for IRAS F01391+51 19 



dilution effect expected in the IRAS data, we conclude that 
the 50% flux discrepancy is due to faint emission in the ex- 
tended nebulosity that was not previously detected by IRAS. 

3.2. Possible Two-Component Structure in the Central Torus 

In this section we give a more detailed analysis of the two- 
peak structure of the central torus. The measured peak sepa- 
ration is 36."7, 58."5, and 42."3 respectively in the 24, 70, and 
160/im bands. Figure|3t displays surface brightness profiles 
along the equatorial plane taken over the whole width (40") 
of the rectangular optical core. 

If the equatorial density profile follows a power law with the 
maximum density at the inner radius (i.e., p = pinf'", where 
Pin is the density at the inner radius, ^n), the largest column 
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Fig. 3. — Surface brightness profiles of the emission peaks along the equa- 
torial plane: (a) profiles at the resolution of each band (FWHM of 6", 18", 
and 40" for the 24, 70, and 160/im bands, respectively) and (b) (convolved) 
profiles at the resolution of the 160/im band. The cuts are made along lines 
of 40" width perpendicular to the bipolar axis (of PA = 45°). Positive offsets 
corresponds to the southwest side of the nebula. While the 70 and 160/xm 
peaks arise from the same locations, the 24/im peaks arise from spatially 
distinct locations closer to the center 



density of the far-IR emitting material along the line of sight 
occurs near the inner radius (at locations ri,, away from the 
center along the equatorial plane). This means that the far-IR 
peaks occur near the inner radius when the dusty circumstellar 
shell is optically thin to far-IR radiation, as has been shown in 
dust radiat ive transfe r models for proto-P N dust shells (e^ , 
Ueta et alJl2 00 1 ; Me ixner et alJl20()2t fUeta & MeixneilBool 
Meixner et al..-2004) . 

When the spatial resolution of the imaging bands is taken 
into account, the peaks are expected to occur progressively 
closer in the bands at longer wavelengths until the peaks be- 
come unresolved. This is why the IQfim peak separation is 
larger than the 160/im peak separation. In fact, when the 
70/.tm image is convolved to have the spatial resolution of the 
160/im band, the 70 and 160/im peaks occur at identical loca- 
tions as demonstrated in Figure l^^- Given the spatial coinci- 
dence of the emission peaks in the 70 and 160/im bands, we 
interpret that these peaks represent the limb-brightened edges 
of a dusty torus. In this picture, the peak separation at the 
70 and 160/im bands corresponds to the inner diameter of the 
dusty torus at a relatively low temperature. 

On the contrary, the 24/im peak separation does not follow 
this trend. If the peaks in the 24/im band represented the Wien 
head of the spectral energy distribution (SED) of the dusty 
torus seen in the 70 and 160/im bands, the peaks should have 
been well resolved with the peak separation slightly larger 
than that in the 70/im band. However, this is not the case. 
Similarly, the 70 and 160/im peaks cannot be the Rayleigh- 
Jeans tail of the SED of a dusty torus whose presence can 



be inferred from the 24/im map. Since the beam size of the 
160/im band is 40", the 24/im peaks (36 "7 separation) cannot 
be resolved. The 24/im band profile convolved with the 40" 
beam indeed shows only one peak as demonstrated in Figure 
ISj). Hence, while the 24/im two-peak structure also suggests 
the presence of an edge-on torus, this torus must possess dis- 
tinct physical properties with respect to the torus seen in the 
70 and 160/im bands. 

Imaging in various optical emission lines have suggested 
that ionization fronts are found only in very selected regions 
in the core and that the lobe edges are not ionization fronts 
(Hahck 1987; Balick et al. 1992). In particular, BaHck ( 19811 
presented an [N II] image showing two emission peaks about 
70" apart in the core and a He II image displaying two peaks 
about 50" apart. The [N II] and He II peaks respectively rep- 
resent "low" and "high" ionization regions. The [N II] peak 
separation is about 45" if convolved with the 40" beam, and 
thus, these peaks spatially corresponds to the possible low- 
temperature dusty torus revealed by the 70 and 160/im maps. 
Similarly, the He II peak separation is about 40" if convolved 
with the 6" beam of the MIPS 24/im band, and thus, the He II 
peaks spatially corresponds to the other torus revealed by the 
24/im maps. This is also demonstrated by the lAC optical 
images in Figure ^ (and Figure |6l for a close-up of the equa- 
torial region). The "high" ionization region represented by 
[O III] emission corresponds to the 24/im peaks, while the 
"low" ionization region delineated by Ha and [N II] emission 
is coincident with the 70 and 160/im emission distribution. 

In fact, Leene & Pottasch ( 1987) have shown that the mid- 
and far-IR structures of the Helix nebula (NGC 7293) are 
markedly distinct in the IRAS data, in which emission in 
the mid-IR (< 25/xm) is concentrated in the central region. 
These authors have attributed the 25/im band emission mostly 
to the [O IV] line at 25.9/im. The following investigations 
have shown that the [O IV] distribution is also co-spatial with 
the H e II distribution in the nebula (O'Dell 1998; Sceck et aO 
Moreover, H2 emission, which suggests the presence 
of a high density molecular reservoir, arises from the limb of 
the bipolar waist , which is spatially coincident with the 70 and 
160/im peaks ( Kastn er et al.ll996tlHora et al.P 2004V 

Therefore, we conclude that the far-IR toroidal structure 
in NGC 650 consists of two distinct components. The first 
component is represented by the 24/im peaks and is a high- 
temperature torus encompassing the regions of high ioniza- 
tion. Given the spatial correlation between the 24/im emission 
region and the emission regions of high ionization species and 
similarities between NGC 650 and the Helix , we also suggest 
that the detected 24/im emission is largely due to the [O IV] 
line at 25.9/im. The second component is represented by the 
70 and 160/im peaks and is a low-temperature dusty torus that 
is exterior to the high-temperature torus and is not yet very 
much ionized by the radiation from the central star. 

3.3. The Remnant AGB Wind Shell 

In this section, we continue our discussion on the low- 
temperature outer torus revealed in the 70 and 160/im bands. 
As reviewed in §1, the circumstellar dust shell of post- AGB 
stars develops its toroidal for m in the innerm ost regio ns dur- 
ing th e superwind phase ('e.g. JMeixner et arfl 999: Uet a et alJ 
1200^. Such a torus would then keep expanding while de- 
creasing its temperature during the subsequent evolutionary 
phases into the PN phase. Hence, we postulate that the dust 
torus seen in the 70 and 160/im maps is the remnant AGB 
wind shell that has been expanding since it was generated by 
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Fig. 4. — Log-log plot of surface brightness profiles along the polar axis 
(top panel) and the equatorial plane (bottom panel). The cuts are made along 
lines of 150" width on images that are convolved with the 40" beam to adjust 
the resolution. The direction of the cuts are indicated at the top corner of the 
panels. The radial offsets are shown from lf'28 to 150" in each direction. 
The slanted line in light gray is a line with the slope of —3 and is given to aid 
readers to read off the slope of the profile. 



the superwind at the end of the AGB phase. 

Given that this outer dust torus is the remnant AGB wind 
shell, its surface brightness profiles would yield valuable in- 
formation pertaining to the AGB mass loss history. If the far- 
thest extent of the highly ionized region has reached only up 
to where the 24/im peaks are, the surface brightness distribu- 
tion beyond the 24/im peaks reflects the mostly pristine AGB 
mass loss history that is still relatively undisturbed by the pas- 
sage of the fast wind and ionization front. 

In Figure 0] we present log-log plots of surface brightness 
profiles along the bipolar axis (top panel) and the equatorial 
plane (bottom panel). These profiles are constructed from cuts 
of 150" width taken in images that are convolved with the 40" 
beam to remove the difference due to resolution. Comparison 
between the top and bottom plots immediately indicates dis- 
tinct mass loss histories along these two directions. 

The equatorial profiles follow a nearly perfect power-law 
with an index of three along the edge of the nebula. In an 
optically thin isothermal dust shell, a power-law density pro- 
file of p{r) oc r~" will yield a power-law surface brightness 
profile of I{r) oc r"*""'*. Thus, the index of three we have 
observed in the equatorial profiles suggests that the equato- 
rial density profile follows r"^. This means that at the end 
of the AGB phase (presumably the most effectively during 
the superwind phase) the rate of mass loss was enhanced 
rapidly, since a constant mass loss at constant velocity yields 
a profile with an index of two. Such a high power-law 
factor in density profiles of the late AGB/superwind shell 
is consistent with recent reports of both obse rvational (e.g., 
iMeixner et al l2()04llHrivnak & Biegingl2()()5l) and theoretical 



(e.g.. lPerinotto et 31.120041 ISchonberner et al.l20'05h studies. 

Comparison between the equatorial profiles (Fig. 0] top) 
suggest possibly different mass loss histories along each edge 
of the torus. While the southwestern profiles remain steep 
with an index of three, the northeastern profiles are less steep 
especially beyond roughly 50" in the 160/im. Thus, there 
appears to be a larger amount of lower temperature dust lo- 
cated at farther away from the central star on the northeast 
side of the torus (i.e., the hmpm/ ho^m ratio becomes large 
as a function of radius along the northeastern edge). One in- 
terpretation is that the period of the increased mass loss was 
less intense yet prolonged on the northeastern edge than on the 
southwestern edge. Such mass loss would result in a higher 
concentration of dust on the southwestern edge of the torus, 
which would cause a stronger dust emission as observed. 

The polar profiles (Fig.|3 bottom) exhibit histories of more 
complex mass loss. The outer part (> 50") appears very sim- 
ilar to the equatorial profiles (showing a power-law with an 
index of three). However, the profiles flatten considerably (to 
an index of less than 0.5) roughly between 50" to 20" and 
then become steep again (to an index of about two) for <, 20". 
Unfortunately, the possible rejevenation of mass loss in the 
inner part is not entirely certain since the inner profiles are 
dominated by the bright torus of 40" thickness. Nevertheless, 
the existence of the plateau region appears real. 

Since the density profile follows in the outer part of 
the shell along both the polar and equatorial directions, it ap- 
pears that the rate of mass loss was enhanced in all direc- 
tions during the AGB phase. From the data, we recognize this 
interesting history of mass loss in this object, in which the 
initially isotropic mass loss enhancement ensued along the 
equator towards the end of the AGB phase while it ceased 
in the polar directions sometime before the end of the AGB 
phase. This is the only apparent change in mass loss during 
the late AGB phase prior to the formation of the equatorial 
density enhancement in the shell. In fact, the observed sur- 
face brightness profiles suggests the presence of a more pro- 
late density distribution around the innermost toroidal density 
distribution. This is consistent with what the layered shell 
model proposes for the structure of the proto-PN shells ('Ueta| 
i2902; Ueta & Meixner 2003) based on the mid-IR and optical 
imagin g surveys of t he proto-PN she lls ( IMeixner et alii 199^ 
lUeta et al...200Q: Meixner et al J2002h as discussed above. 

Therefore, we conclude from these profiles that (1) mass 
loss was enhanced nearly isotropically towards the end of the 
AGB phase, (2) mass loss in the polar directions precipitously 
decreased prior to the end of the AGB phase, and (3) the 
density distribution transformed from prolate spheroidal to 
toroidal during the late AGB phase. Although it is not entirely 
clear which parts of the profiles are pertinent to the superwind 
epoch, it appears evident that the superwind epoch at the end 
of the AGB phase is associated with a cessation of mass loss 
in the polar directions rather than an enhancement of it in the 
equatorial directions. This is a strikingly contrasting view of 
the late AGB mass loss with respect to a commonly-held as- 
sumption that there is an enhancement of mass loss into the 
equatorial plane. Nevertheless, the present data provide an- 
other piece of evidence that supports the presence of a transi- 
tional spheroidal shell between the outer spherical AGB shell 
and the inner toroidal superwind shell as proposed in the lay- 
ered shell model (Ueta 2002). 

A particularly important aspect of the present work is that 
this transitional spheroidal shell has been observed in a bipo- 
lar PN. This is a direct observational evidence for the presence 



Spitzer/MIPS Imaging of NGC 650 



7 



of the spheroidal density distribution beyond the inner torus in 
a bipolar PN. It has been typically (and probably naively) as- 
sumed that there is an absence of matter along the polar axis 
in bipolar shells of evolved stars, proto-PN shells and PNs 
alike. However, the presence of matter along the polar axis is 
evident in this bipolar PN, further increasing our confidence 
in the layered shell model. The present data also demonstrate 
that one can probe the AGB mass loss history via far-IR imag- 
ing of PNs as long as the shell structure is not altered by the 
passage of the ionization front and a fast wind. 

In the following we estimate some quantities concerning 
the remna nt AGB wind shell. By adopt ing the distance of 1.2 
kpc (e.g., iKoornneef & Pottaschlll998l) . the peak separation 
of 42 "3 (at 160/xm) translates into 7.6 x 10'^ cm. Given the 
spatial resolution effect, this value should be taken as a lower 
limit for the inner diameter of the dust torus generated by the 
AGB mass loss. If we adopt a typical expansion velocity of 
20 km s~' (of a superwind) and assume constant-velocity ex- 
pansion, the dynamical age of the dust torus is estimated to be 
6.0 X 10^ yr This value is consistent with theoretical predic- 
tions of the post- AGB evolution into the PN phase of interme- 
diate mass (^ 3 M©) stars (e.g.. Blocker 1995). 

Assuming that there is no significant line contamination in 
the measured fluxes and that the dust shell is optically thin at 
far-IR, we can estimate the dust temperature in the torus by 
fitting the far-IR SED of the torus using the relation. 



/V = ribeamTi/Q ( — ) ^[/(Tdust), 



(1) 



where /v is flux density at frequency i^, fibeam is the instru- 
mental beam area, t{i^)(= T,^a(y/i^o)^) is an optical depth at 
(scaled by a power-law), and Tdust is the (representative) dust 
temperature along the line of sight. Assuming fluxes long- 
ward of 100/im represent the Rayleigh-Jeans tail of the SED, 
the overall dust emissivity index can be fixed to 1.6 ±0.5. 
By setting /3 = 1 .6 in the above relation, we obtain the least- 
squares fits of Tdust and T(160^m) under the two-component 
dust assumption, using the present MIPS fluxes and the past 
IRAS and IRAC fluxes. The dust temperatures are found to 
be 32 ± 14 K and 140 ± 53 K and the optical depth at 160^m 
is found to be 1.2 x 10"^. The SED and fitted curves are dis- 
played in Figure |5] Given the power-law emissivity of the 
dust grains, the optical depth of the shell in the optical (at V) 
is 0.1: the torus is indeed optically thin to visible light. 

If the overall dust opacity, k, is known, the measured far- 
IR flux can be used to estimate the mass of dust in the torus 
through 



(2) 



Here, is the distance to NG C 650, c is the sp eed of light, 
and k is Boltzmann constant. iSemenov et alJ (l2603) have 
computed Rosseland and Planck mean opacity tables for the 
dusty medium based on the dust composition model for ac- 
cretion discs consisting of various silicates, organics (CHON 
material), amorphous water ice, FeS, and iron having vari- 
ous particle types. By adopting k = 0.4 to 0.7 cm^ g"' for 
the C/O = 0.43 and Fe/(Fe + Mg) = 0.3 (i.e., "normal" iron 
abundance) case from their opacity table, the above equation 
yields 1.4 to 2.5 x 10"^ M© for the measured 160/im flux of 
4.83 Jy and the estimated dust temperature of 32 K. With a 
canonical gas-to-dust ratio of 150, the total mass of the shell 
is estimated to be 2 to 4 Mq. The gas-to-dust ratio is highly 
uncertain, and so is the estimate of the total mass of the shell. 




Fig. 5. — The far-IR spectral energy distribution of NGC 650. Fluxes are 
taken from MIPS (square), IRAS (circle), and IRAC (triangle) observations. 
The vertical bars indicate an error in each measurement. The dashed line 
shows the best-fit to the colder AGB component of the dust shell (at 32 K), 
whereas the dotted line shows the higher temperature component (at 140 K). 
The gray solid line indicates the total flux. 



Nevertheless, the estimate is consistent with the inferred ini- 
tial mass of the central star (^ 3 Mq) based on the dynamical 
age of the torus compared with theoretical post-AGB tracks 
(see above). Given the extent of the observed far-IR emission 
(~ 50") and the inner radius of the torus (^ 20-30") plus the 
assumed wind velocity of 20 km s"' , the above values yield 
an overall mass loss rate of 2-5 x 10"^ M© yr"'. 

3.4. The Post-AGB Shell 

We will now focus our attention to the inner, high- 
temperature torus seen in the 24/im band. The 24^m peaks 
do not simply represent the inside wall of the dust torus 
which harbor high-temperature matter that is illuminated by 
the central star. If this were the case, the 24/im peaks should 
have been spatially coincident with the 70 and 160/im peaks. 
Based on a spatial correlation between the 70 and 160/im 
peaks and the high ionization regions of He II and [O III], we 
have concluded that the 24/im emission is probably largely 
affected by the [O IV] Hne at 25.9/im. 

It is not impossible, however, to have dust emission in such 
highly ionized regions. For example, the past Infrared Space 
Observatory observations detected rather strong dust emission 
at 90/im in the central region (i.e., the most ionized region) of 
the Helix nebula ( Speck et al. 2002). The 70/im map indicates 
that there is another emission peak near the center of the neb- 
ula (Figure|3t; right panel of Figure|6j. Even in the 160/im 
map, there appears such emission at the nebula center: the 
"trough" between the emission peaks is very shallow (bottom 
panel of Figure^ Figure|3}. 

Aperture photometry of this central region (within the in- 
ner radius of the low-temperature torus) yields 0.54 and 0.035 
Jy respectively in the 70 and 160/im bands. These values are 
surprisingly close to the fluxes of the high-temperature SED 
component at these wavelengths (0.51 and 0.031 Jy respec- 
tively at 70 and 160/im; Figure |5}, which cannot at all be 
entirely attributed to photospheric emission of the central star 
(3 X 10"'' and 6 x 10 Jy respectively at 70 and 160/im extrap- 
olated from t he photospheric photometry in the IRAC bands 
presented by Hora et al. 2004). Using the estimated far-IR 
fluxes from the inner torus and k ^ 0.2 cm^ g~' for dust of 
140 K from the opacity table of ISemenov et al J ( [20031) in eqn 
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Fig. 6. — The deconvolved 24 (left) and lOfim (right) MIPS maps of the central 100" X 100" region (N is up. E to the left) in contours overlaid with the lAC 
optical maps in the [O III] (left) and Hq+[N II] (right) bands in grayscale. The contours are, from top to bottom, 90 to 10% of the peak with the 20% interval. 
The location of the central star ([0,0] offset) is indicated by a "star", showing the location of the lOfim central peak with respect to that of the central star. The 
other display convention follows that of Fig. HI 



(2), the mass of dust in this central peak is 2.5 X 10"^ Mq. 
Given that only part of the measured 24/im flux is dust-origin, 
this mass estimate can only be a rough upper limit. The inner 
torus is thus very tenuous. 

Figure |6l shows the central 100" x 100" region of the de- 
convolved 24 and 70/im contours overlaid respectively with 
the [O III] and Hq!-i-[N II] grayscale maps. The deconvolved 
24/im map shows an extension in the southwest peak towards 
the center of the nebula. This extension indicates the presence 
of some matter that fills the inner cavity of the torus. The de- 
convolved lOfim map also indicates the presence of some dust 
in the inner cavity of the torus, which probably represents the 
Rayleigh- Jeans tail of the high-temperature dust component 
as has been discussed above. 

Thus, the inner cavity of the torus appears to be filled with 
highly ionized gas and some surviving dust behind the ion- 
ization front that has already engulfed the inner torus. Where 
does this matter in the inner cavity come from? There are two 
possibilities: this matter could come either from the central 
star or back from the torus. If the cavity-filling matter comes 
from the star, the matter is likely part of a fast wind. Then, 
it is very peculiar that a fast wind has generated wind-blown 
lobes along the bipolar axis, while maintaining a rather high 
concentration of matter along the equatorial plane. 

Alternatively, matter in the torus that is engulfed by the ion- 
ization front can flow back towards the central star and fill the 
inner cavity of the torus, if the pressure from a fast wind can- 
not hold it. The back-flowing material from the torus initially 
fills the cavity, but eventually disperses into the lobes. The 
cavity-filling matter can get ionized by the radiation from the 
central star to generate even 0+++. In fact, (iBrvce et all 19961) 



have observed multi-component [O III] 5007 emission line 
profile in the central region of the nebula, and the low-velocity 
components can be interpreted as back-flowing ionized gas. 
Therefore, the inner torus represented by the 24/im emission 
seems to consist of the ionized portion of the torus behind the 
ionization front, some ionized matter that has flown back to- 
wards the star filling the inner cavity of the torus, and some 
surviving dust species. Thus, the presence of the post-AGB 
wind shell in NGC 650 is not very likely. 

4. CONCLUSIONS 

We have investigated the AGB and post-AGB mass loss his- 
tory in a bipolar PN, NGC 650, using far-IR maps at 24, 70, 
and 160/im taken with MIPS on-board the Spitzer Space Tele- 
scope. The far-IR surface brightness distribution shows two 
emission peaks at all three MIPS bands. These peaks repre- 
sent the limb-brightened edges of a near edge-on, optically 
thin dusty torus. Hence, this PN's far-IR structure is con- 
sistent with the known mid-IR structure of the PN progeni- 
tors, suggesting that a dusty torus in the circumstellar shells of 
evolved stars endure the transition from the post-AGB phase 
to the PN phase. The measured fluxes above three cr are 
4.51 ±0.04, 6.04 ±0.30, and 4.83 ±0.95 Jy respectively at 
24, 70, and 160/im. 

While the 70 and 160/im peaks are spatially coincident, the 
24/tm peaks are found interior to the peaks in the other bands. 
Based on the spatial correlation with the distribution of [N II], 
Ha-i-[N II], He II, and [O III] in the optical, we have concluded 
that the far-IR toroidal structure consists of two components. 
The 70 and 160/tm peaks represent low-temperature 30 K) 
dust continuum emission from the remnant AGB wind shell. 
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The 24/im peaks correspond to the [O IV] Hne emission at 
25.9/im arising from the highly ionized part of the inner torus 
which has been engulfed by the ionization front and the ion- 
ized matter that has flown back from the torus and fiUed the 
inner cavity of the torus. 

The dynamical age of the remnant AGB wind shell is 6 x 
10^ yr, assuming a constant expansion velocity of 20 km s"' 
and the distance of 1 .2 kpc. The least-squares fits to the far-IR 
SED have yielded the temperatures of 32 ± 14 and 140 ± 53 K 
for the two-component dust shell system and the optical depth 
of Tieo^m = 1.2 X 10"^ (i.e., Ty = 0.1). The dust masses of the 
remnant AGB wind shell and the inner torus are estimated to 
be 1.4 to 2.5 x 10"^ and 2.5 x 10"^ Mq, respectively. For a 
canonical gas-to-dust ratio of 150, the total mass of the shell 
is 2 to 4 M0, yielding an overaU mass loss rate of 2-5 x 10"^ 
Mq yr-1. 

From the far-IR surface brightness profiles of the nebula 
that show a power-law index of as steep as three, we have also 
reconstructed this PN's AGB mass loss history in which (1) 
mass loss was enhanced nearly isotropically towards the end 
of the AGB phase, (2) the enhanced mass loss in the equato- 
rial directions ensued at the end of the AGB phase while that 
in the polar directions precipitously ceased prior to the end of 
the AGB phase, and (3) the density distribution transformed 
from prolate spheroidal to toroidal during the late AGB phase. 
Hence, the superwind phase at the end of the AGB phase is as- 
sociated with a cessation of isotropically enhanced mass loss 
in the polar directions rather than an enhancement of it only 
in the equatorial directions. 

This is a strikingly contrasting view of the late AGB mass 
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loss with respect to a commonly-held assumption that there is 
an enhancement of mass loss only along the equatorial plane. 
The present data thus provide another evidence for the pres- 
ence of a transitional spheroidal shell between the outer spher- 
ical AGB shell and the inner toroidal superwind shell, further 
increasing our confidence in the layered shell model for the 
circumstellar shells of post- AGB stars. Moreover, the data in- 
dicate that the mass loss enhancement into the equatorial di- 
rection is not strictly uniform in all azimuthal angles. Mass 
loss appears to be less intense but more prolonged on the 
northeast edge of the torus than the southwest edge, yield- 
ing higher concentration of dust at the inner radius on the 
southwest edge. Nevertheless, the present data also demon- 
strate that one can probe the AGB mass loss history via far-IR 
imaging of PNs as long as the shell structure is not altered by 
the passage of the ionization front and a fast wind. 
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